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Helicity-dependent incoherent pion photoproduction in the reaction 7c! — > nNN is studied in the 
framework of the diagrammatic approach. Contributions to the reaction amplitude from diagrams 
corresponding to impulse approximation as well as AW and nN interactions in the final state 
have been evaluated. The elementary *yN — > nN operator is taken from the MAID and SAID 
models. A detailed comparison of the predictions with recent experimental data by the GDH and 
A2 collaborations at energies below 500 MeV is presented. Reasonable agreement with the data 
on the yields and cross sections for polarized beam and polarized target has been achieved in all 
channels. The unpolarized data of the GDH and A2 collaborations have also been analyzed within 
the approach. A strong overestimation for the neutral channel has been found. At the same time, the 
model provides a quite satisfactory description of the unpolarized data for the charged channels. The 
sensitivity of the obtained results to the choice of the elementary "/N — Y ttN operator is discussed 
in detail. The contribution of the 'yd —¥ ttNN reaction to the GDH sum rule for the deuteron up to 
a photon energy of 1.65 GeV has been evaluated with the result of 235 ± 25 yh. 
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I. INTRODUCTION 

For more than a decade, the GDH and A2 collaborations have been carrying out intense experimental studies aimed 
at the verification of the famous Gerasimov-Drell-Hearn (GDH) sum rule A comprehensive overview concerning 

the status of these experiments is given in Ref. 0]. 

The GDH sum rule gives an integral relation between anomalous magnetic moment k of a spin-5 particle and the 
difference of the total photoabsorption cross sections with parallel, op, and antiparallel, a a, photon and target spin 
alignments and reads 
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where E n is the photon laboratory energy, m is the particle mass, and a = 1/137. Equation (jlj shows that ground- 
state properties given by n and m are related to the energy-weighted excitation spectrum of the particle. The GDH 
sum rule relies on the basic physics principles of Lorenz and gauge invariance, unitarity, crossing symmetry, and 
an unsubtracted dispersion relation applied to the forward Compton amplitude. Therefore, measurements of the 
right-hand side (rhs) of this equation can serve as a fundamental cross-check of these principles. 

The left-hand side (lhs) of Eq. (jJTJ) for the proton is 205 /ib. The integrand in the rhs of Eq. ([TJ at energies from 
0.2 to 2.9 GeV was determined by the GDH and A2 collaborations in a series of experiments with circularly polarized 
photons and longitudinally polarized protons 0-0- Involving theoretical predictions from threshold to 0.2 GeV and 
above 2.9 GeV, the GDH group obtained for the rhs of Eq. (JTJ) 212 ± 6(stat) ± 16(syst) fib, which is in good agreement 
; "*] ■ with the GDH sum rule prediction. 

The check of the isospin structure of the GDH sum rule requires measurements on the neutron. However, because 
of the absence of a stable, dense, free-neutron target, direct measurements for the neutron are impossible so that 
information on the neutron can be extracted from deuteron data. Experiments with circularlypolarized photons 
and longitudinally polarized deuterons were also performed by the GDH and A2 collaborations |7Hllj| as well as by 



the LEGS Collaboration [12j. In fact, however, a contribution to the neutron GDH integral from deuteron data was 
extracted only in Ref. Q and therewith in the limited energy range of 815-1825 MeV. The extraction procedure relied 
mainly on an assumption that at these energies the incoherent, quasi-free meson production reactions dominate; i.e. if 
one supposes that the equality A<7d = Aa p + Aer„ (Act = ap — a a) is valid. Such an assumption seems to be in general 
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agreement with the results obtained by Arenhovel, Fix, and Schwamb (AFS) [ljj. Using theoretical predictions for 
the unmeasured energy regions, the authors of Ref. obtained the value of 226 /A> for the rhs of Eq. (fTJ) which is 
in agreement with the GDH sum rule value of 233 fib. Of course, this agreement should not be considered as an 
experimental verification of the GDH sum rule for the neutron, because the net experimental contribution to this 
value is only about 15%. 

The helicity-dependent cross section on the deuteron in the energy range 200 < E 1 < 800 MeV was also measured 
in Ref. [lol ] . Combining the obtained data with those from Ref. @ and using again the assumption on the incoherency 
of the proton and neutron contributions in the range 200 < E~ < 1800 MeV, the authors of Ref. [T3| obtained the 
value of 197 fib in this range. The said assumption seems to be rather rough below, say 500 MeV, where the 7c? — > np 
and the yd — > w°d reactions contribute to the polarized total cross section, nevertheless, the extracted value indicates 
that the GDH integrals for the proton and neutron are of the same order of magnitude. 

Together with the verification of the GDH sum rule for the neutron, the helicity-dependent cross sections on the 
deuteron are also required to test this sum rule for the deuteron itself. The anomalous magnetic moment of the 
deuteron is very small (K4 — —0.143 nm) because, first, the anomalous magnetic moments of the proton and neutron 
are almost equal in magnitude but opposite in sign and, second, nucleon spins in the deuteron are predominantly 
parallel. In view of this, the lhs of Eq. (fTJ for the deuteron is 0.65 /ib, which is more than two orders smaller than the 
nucleon values. Therefore, there should be large negative contributions almost totally canceling the nucleon contribu- 
tion. This point was discussed for the first time in Ref. 14] . It was found that a huge negative contribution amounting 
to about —400 fib indeed exists and stems from the deuteron photodisintegration reaction. Later on, a comprehensive 
analysis of different channels in the deuteron GDH integral was performed within the AFS framework [13j]. They 
evaluated effects from deuteron photodisintegration and single and double-pion production as well as from single-?? 
production with the total result of 27.31 fib. 

If one intends to extract the neutron value for the GDH integral from deuteron data or verify the deuteron GDH 
sum rule, in either case one has to be confident that theoretical models used in analyses of these data provide reliable 
descriptions of the helicity-dependent cross sections for different reactions on the deuteron. As was found in Refs. 
fl5j |. a large contribution to the deuteron GDH integral stems from incoherent single-pion production 7c? — > nNN. 

Intense experimental investigations of the reaction yd — > ttNN have recently been performed by the GDH and A2 
collaborations (9l4TTI| . References [§, [l(| have provided data on the helicity-dependent total inclusive and semiexclusive 
cross sections in the energy range from 200 to 800 MeV. In the very recent work of the same collaborations [ll[ , the 
helicity dependence of the differential cross section in the yd — > irNN channel has been measured for the first 
time in the A region. There are also measurements of the helicity-dependent total cross section of the reaction 
yd—> ir°X 0, [Hj], where X — np (the incoherent channel) or X — d (the coherent channel). 

The main goal of the present work is to analyze the data of Refs. 0-[IH on the yd — > ttNN reaction in the framework 
of the diagrammatic model built in Ref. [16]. In that work, we restricted ourselves to the discussion of unpolarized 
and single spin-dependent observables due to the lack of data on double spin-dependent observables at the time when 
the work was being prepared. It was found that the model [16| provided a quite satisfactory description of all available 
data except for the unpolarized differential and total cross sections in the yd — > ir°np channel for which the data were 
notably overestimated. In the present work, we compare predictions of the approach [16] to the polarized data from 
Refs. @4lll w ith special emphasis on the model dependence of our results. Sources for theoretical uncertainties are 
discussed in Ref. [Tgl, to which the interested reader is referred. 

Another principal goal of this work is to evaluate the contribution from the yd — > ttNN channel to the deuteron 
GDH integral. This task requires an extension of the model up to about 1.5-2 GeV. Although in Ref. [l6| we 
restricted our consideration to the A resonance region, the model is applicable at higher energies, too. As was shown 
in Ref. [13] , one can expect that in the framework of the diagrammatic approach, contributions due to NN and ir N 
final state interaction (FSI) are reliably calculated up to a photon energy of 800 MeV. Taking into account that these 
contributions are small above 500 MeV, they can obviously be disregarded for E 1 > 800 MeV. Hence, the upper limit 
for the applicability of the model depends on that for the impulse approximation (IA) which, in turn, depends on 
the applicability domain of an employed elementary operator yN — > irN of pion production on the nucleon. In the 
framework of the approach [l6| , this operator is taken from multipole analyses SAID [l8l420| and MAID [2l| and can 
be used up to 2 GeV for the first analysis and up to 1.65 GeV for the second one. 

The contribution from the yd — > ttNN channel to the GDH integral up to E 1 = 1.5 GeV was evaluated in the AFS 
approach [131 ] . but the sensitivity of obtained results to the yN — > itN operator was not investigated in that work. 
An analysis of Ref. [l5| shows that the value of the integral at integration up to 350 MeV manifests its considerable 
dependence on the elementary operator. Because the energy domain from threshold to 350 MeV provides only about 
a half of the total contribution of this reaction to the deuteron GDH integral, it is of high interest to recognize how 
sensitive the total integral is to the yN — > ttN operator. In the present paper, we evaluate the integral up to a photon 
energy of E 1 = 1.65 GeV and investigate a question of possible theoretical uncertainties for its value. 
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The work is organized as follows. In Sec. UH the kinematic relations used in the calculations as well as definitions 
for observables are reviewed. A brief description of the theoretical model and its ingredients is given in Sec. Mil 
Section HVl contains results for the unpolarized and helicity-dependent yields and cross sections and their comparison 
with the data available in the considered kinematic region. Results of the evaluation of the jd — > irNN channel 
contribution to the deuteron GDH integral are also given in Sec. IIVI In Sec. [V] we summarize the main conclusions 
and results. 



II. KINEMATICS 



Actual calculations of the reaction amplitude and observables are done in the laboratory frame where k = 
(E 7 ,k), pd — (M, 0), Ptt — (e^q), pi — (et,pi), and P2 = (£2,1*2) are the four-momenta of the initial photon 
and deuteron and the final pion and nucleons, respectively. For the case when the final pion is detected, we take as 
independent kinematic variables the photon energy E 11 the value of the pion momentum q = \q\, the pion polar angle 
Gtt, and the solid angle fig (05 , <t>s ) of the relative momentum Pmn of the final nucleon-nucleon pair. These 
variables totally determine the kinematics. 

The unpolarized differential cross section is given by 

mr='™ I £ k«i^>i 2 , (2) 

where {m2m\ \T\Xmd) is the reaction amplitude, and the phase space factor is 

fNN ~ J^f 4£ 7 e„WW (3) 



with Wnn — \J (pi + P2) 2 and m is the nucleon mass. Spin states of the two nucleons and deuteron are 7712, mi, and 
md, respectively, and they are chosen with respect to the z axis, which is defined by the photon momentum k. The 
symbol A stands for the photon helicity. An extra factor of 1/2 must be included in the rhs of Eq. (J3J in the case of 
charged pion production. 

We will be also interested in double polarized observables for the polarized photon beam and polarized deuteron 
target. These are the parallel, P, and antiparallel, A, cross sections defined as 

J^- = fNN \ ^ (K^im + l + ^P + K^milTl-l-l)! 2 ), (4) 

- Pnn ni2mi 



-Pa. 

da A 



dqdttg 



= fNN ~J1 (k^mjTl + l-l)! 2 + |(m 2 m 1 |T|-l + l)| 2 ). (5) 



To obtain the semi- inclusive differential cross sections da/dfl^, dap/dfl^, and daA/dfl^, the right-hand sides of 
Eqs. and ([5]) have to be integrated over q and the solid angle Qp NN , i-e., 

da f 9 , f da , . 

dn„ = J q ,J NNqdq J dn ^dqm—^ (6) 

! ' fNN q 2 dq I dttp MM -j^^ — , (7) 

'nn 



dn n Ja^in J ^ NN dqdflp 



^ f^N 1 2 dq ( dttp — . (8) 



dfl^ J ami n J *nn dqdttn 

The upper and lower integration limits in the laboratory frame arc defined by the following relations 

.max _ „max (9i) = 1 \ £ + ^ + M) ^2 _ bfJ 2 1 ) (9) 



q =q 



q mm = g mm 



(9^) = maxjo, - aE^z„ - (E y + M)y^a 2 - bfj? }, (10) 



where a = (W 2 d — 4m 2 + ^ 2 )/2, b = (E 1 + M) 2 ~ E 2 z 2 , = cos 9^, and ^ is the pion mass. 
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For the neutral channel, the emitted proton was registered in the experiment [Tl(. In this case we chose as 
independent kinematic variables the photon energy E 11 the value of the proton momentum p p = \p p \, the proton 
polar angle P , and the solid angle Hp (Op , <fip ) of the relative momentum P w „ of the final meson-neutron pair. 
Equations analogous to Eqs. (HJ and (fTO)) are now of the form 



.max max 



^ laX (©p) 

p; nin = pf n (e P ) 



1 

B 

1 



AE 7 z p + [E 1 + M)V 'A 2 - Bra 
{jo, AE lZp - (E 1 + M)\J A 2 - Bra 2 } 



(11) 
(12) 



where A = (W 2 d - 2mfi - ^. 2 )/2, B = (E 7 + M) 2 - E 2 z 2 , and z p = cos6 p . 

Only a fraction of the total interval [p p mn ,j> p aax ] was accessed in conditions of the experiment [ill ]. There were 
registered the protons with the absolute value of the momentum p p > , where 



fCp - PS?(Q,) = [210 + 0.01 x (6 P - 90°) 2 ] MeV/c, 



(13) 



so that the differential 7T°np yields Y, as they were called in Ref. [TT|, were measured rather than the differential cross 
sections. In analogy to Eqs. ©-(jH) they are defined as 



dY 

dY P 
dQ p 

dY A 
dQ n 



J Un P 2 P dpp J dflp^^ i ^ \(m p m n \T\\m d )\ 2 , 

j Un P 2 P dpp j dflp^ n i ^ (\( m p m n\T\ + 1 + 1)| 2 + \(ra p ra„\T\ -1-1) 

Pcxp m p Tn n 

j Un P 2 P dpp j dQp^^ ^ (\(m p m n \T\ + 1 - 1)| 2 + \(m p m n \T\ - 1 + 1) 



where 
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(14) 
(15) 
(16) 

(17) 



and W.„ 



\APtt +Pn) 2 



III. THE FORMALISM 



In this section we briefly outline the formalism used to evaluate the jd — > ttNN reaction amplitude. It was described 
in Refs. [16l . |22j, |23j, to which the reader is referred for more details. The diagrammatic approach is exploited to 
calculate the 7c? — > irNN reaction amplitude. Because in this work we are not interested in the threshold region, we 
reduce the set of diagrams in comparison to that considered in Ref. [23|. Specifically, only contributions from the 
three diagrams shown in Figs. [IJa)-(c) to the reaction amplitude have been taken into account. Diagram[IJa) is often 
refereed to as IA. Sometimes it is called the pole diagram. Apart from it, the diagrams [TJb) and (c) corresponding, 
respectively, to NN and irN rescattering in the final state have been considered in the present calculation. 




FIG. 1: Diagrammatic representation of the 'yd — > tvNN amplitude. Diagram (a) corresponds to IA. Mechanisms with NN and 
irN rescattering in the final state are shown in diagrams (b) and (c), respectively. Diagrams with the permutation Ai A2 
are not shown. The spectator nucleon indicated by an x is on-shell. 
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One remark should be made here. A detailed analysis of the helicity dependence of the —> ttNN reaction up the 
A(1232) resonance is presented in Ref. [lfij]. Special attention was given in that work to the near-threshold region. 
Within the approach [15| the same set of diagrams as presented in Fig. [1] was taken into account. One can expect 
such a model to be quite realistic near threshold for the descri ptio n of the charged channels. But this is not obviously 
to be the case for the neutral channel. As was shown in Ref. [23|], a two- loop diagram as in Fig. Ufa) which includes 
simultaneously TtN rescattering in the intermediate state and np rescattering in the final state has to be taken into 
account at threshold energies. A similar mechanism is known to be very important near threshold for coherent ir° 
photoproduction on the deuteron — > TT°d [24-26j . 




FIG. 2: Two-loop diagrams of the reaction 'yd —¥ nNN. The permutation iVi o N2 is not shown. 



The importance of diagram (a) in Fig. [5] at threshold energies is that it contains, together with the np scattering 
amplitude, a block with charged pion photoproduction from the nucleon. The dominating threshold electric dipole 
amplitude Eq + for the charged channels is about 20 times larger in absolute numbers than that for the neutral channels. 
For analogous reasons we also expect that a two-loop diagram Hfb) with NN rescattering in the intermediate state 
and ttN rescattering in the final state may also be important near threshold. Therefore, the fact that the model 
including only three contributions shown in Figs. [IJa)-(c) can be used for the description of the jd — >• ir°np reaction 
at threshold energies is not obvious, and further study is needed to clarify this point. There is one more two-loop 
diagram, shown in Fig. [He), which includes two blocks with irN scattering. Because the irN scattering lengths are 
about two orders smaller than those for NN scattering, the role of this diagram in the threshold region is expected 
to be much smaller than that of diagrams [2]^ a) and [2jb) . 

One can anticipate the effect from the above two-loop diagrams to decrease together with increasing E 7 when 
the mentioned enhancement due to the charge pion exchange disappears. Nevertheless, to be fully confident in the 
smallncss of such contributions, their explicit evaluations are needed. This task, however, is expected to be very 
difficult, because as we have learned evaluating diagram [2Ja) in Ref. [23|], practical calculations are extremely time 
consuming especially when the inclusive processes are considered, and, therefore, integrations over momenta of final 
particles have to be carried out. Note that in the case of the exclusive channels, the kinematic regions exist where the 
diagrams in Figs. [TJa)-[lJc) can be suppressed, and other mechanisms, in particular, the two- loop diagrams [2fa)-[2jc) 
become important and have to be taken into account (see Refs. [27:, 28] for detailed discussions of this point). 

Let us now write explicit expressions for the amplitudes corresponding to diagrams [T](a)-{TIc). One has for the IA 
amplitude 

(m 2 mi|T IA (p 2 ,Pi,g;fc)|Am d ) = ^ *™^ m (p 2 ) (mi|T 7 ;v_ MrJ vi(pi, q; A;)|Am), (18) 

m 

where ^^im ^ s ^ ne deuteron wave function having the form 

*ZU(P) = E i L C 1 &^Cfr B ^ L Yr>&u L <p). (19) 

L=0,2 

In Eq. (fH)]) . Y™ L (p) are the spherical harmonics, and Cj 1 ^ j 2 m 2 are ^ ne Clebsch-Gordan coefficients. The S and 
D wave function amplitudes uq(p) and 112 (p) are taken for the CD-Bonn potential (29[. One should emphasize that 
calculations with another realistic version of the nucleon- nucleon potential, namely, with the Nijmegen model [30] , have 
led to essentially the same results, so we have found no sensitivity of our predictions to the choice of a potential. Note 
that in Eq. (1181) and in the amplitudes that follow, we do not write explicitly those corresponding to the permutation 
Ni «-» ./V2, but they are included in the calculations. 

The symbol T^n^^n in Eq. (118[) stands for the elementary operator of pion photoproduction on the nucleon. As 
is explained in detail in Ref. [16[ (see, also Ref. [ll[), ambiguities in predictions for observables in the 7c? — > -kNN 
reaction stem mainly from the manner in which T^n^ttN is embedded into the deuteron. In all recent calculations II 5l— 
[ItJ , this operator is taken in the on-shell form and it is parameterized using either the multipole analyses SAID |l8| 
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and MAID or the effective Lagrangian approach (ELA) (3l|. Although the models are close in their predictions 
for observables in the reaction 7^ — > 7riV, they are not identical. It is evident that this difference will manifest itself 
in results of the evaluation of Eq. (TT8l) . To estimate ambiguities caused by different choices of the on-shell operator 
T 7 jv->ttJV) we perform our calculations with the latter taken from the recent multipole analyses SAID [l8j (solution 
SP09K) and MAID07 (H. 

Furthermore, the on-shell operator T^n^^n depends on four invariant amplitudes [33[ for which there exist different 
options that are equivalent in the on-shell case. This equivalence is broken when one or both nucleons are off their 
mass shells. Just such a situation takes place in Eq. (TT5)) where the nucleon TV is off its mass shell. In this work, we 
evaluate the operator T^jy-y^N with two sets of the invariant amplitudes, Ai and as they are called in Ref. and 
where their definition can be found. At least, this allows one to estimate the possible uncertainties in the predicted 
observables due to off-shell effects in T 7 jv->7riv- 

The matrix element corresponding to Fig. [ljb) reads 

(m.rmlT^ip^q-MXma) = mj g*', ^ _ ^ + , Q 

x ^ (po U t,m 2 m 1 \T NN (E)\p in ,m s m')(m s m / \T 1A (p s ,p',q; k)\Xm d ). (20) 

m s m' 

Here T IA is the IA amplitude given by Eq. (| . The relative momenta of the N2N1 pair after and before scattering are 
Pout = (P2 ~ Pi)/2 and pi n — p s — (f>2 + pi)/2, respectively, and E = p^ ut /m. The half-off-shell NN scattering matrix 
Tjv n has been obtained by solving an equation of the Lippmann-Schwinger type for the CD-Bonn potential [29[ . All 
states with the total angular momentum J < 2 have been retained in T/vat. 

In the evaluation of Eq. (|20p there is another ambiguity due to different prescriptions for the determination of 
the total invariant energy entering the 7 AT vertex. As was proposed in Refs. [U HH, we define the latter as if the 
spectator nucleon N s was on its mass-shell. Other prescriptions can also be found in the literature. For instance, in 
Ref. [13] the active nucleon N is put to be on its mass-shell. We have evaluated Eq. (|2"0")) with these two prescriptions 
and found only at most a 2% variation for all observables considered in the next section. This means that the above 
ambiguity can be safely disregarded. 

The matrix element corresponding to Fig. [He) with irN rescattering in the final state is 

X 5Z [(P2i m 2\ T nN( W *N 2 )\Ps,rn s )(m s mi\Tlo Nl {Ps,P^ (21) 

rn s 

~ (P2, m 2\T^( W 7vN 2 )\Ps, m s)(m s m 1 \Tlt Nl (Ps,PiA^i;k)\Xm d ) , 

where the asterisk denotes variables in the KN2 center-of-mass (cm.) frame, i.e. p* = — q*, where p* and q* are 
the cm. initial nucleon and pion momenta, respectively, in the -kN scattering vertex. Analogously, p 2 * = —q*, 
where p 2 * an d q* are the cm. final nucleon and pion momenta, respectively. The momenta with the asterisk and 
without it are related to each other through a boost transformation with the velocity {q + P2)/ + £2)- We take 
the spectator nucleon N s to be on its mass shell, i.e. e* = \Jp*^ + m 2 . The total energy W^n 2 of the ttNz pairs is 

W* Na = y/(q + P2) 2 = e% + 4 = + m 2 + V<1* 2 + M 2 - 

The half-off-shell irN scattering matrix T^n has been obtained by solving the Lippmann-Schwinger equation for a 
separable energy- dependent ttN potential built in Ref. 1341 . We have checked that the results remain essentially the 
same as the T^n found from a meson-exchange model [3a. l36j |. The factor i^^vd^rJ) m Eq. (EU is introduced to 
take into account the off-shell nature of the intermediate meson 71^ . In accordance with Ref. [3J] , we choose the form 
factor FnN(q) of the ttNN vertex to be of the monopole form 

f **(q) = "nr~ 2> ( 22 ) 

K + T 

with a cutoff parameter A^ = 650 MeV/c. 

The superscript indices "0" and "ch" in the T^n amplitude in Eq. (f2Tj) stand for the neutral, 7r°A^ — > tt°Ni, or 
charge exchange, 7r ch iV — > 7r A^i, channels, respectively. The IA amplitudes for tt° photoproduction, T*$ N , and for 

charged pion photoproduction, T^ hN , are given by Eq. (fTHjl . 
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IV. RESULTS AND DISCUSSION 



A. Unpolarized yield and cross section 



We begin the discission of the results with the unpolarized yield and cross section. Special emphasize is given 
to their comparison with the very recent experimental data from Refs. [ij-[TTl|. A detailed comparison with older 
data can be found in Ref. [16| . As a rule, we will not compare ourpredictions with the results of other theoretical 
approaches [H, E3) H3| because this has already been done in Ref. [ll| . 



80 r- 
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FIG. 3: Unpolarized differential yield for the -fd — > ir°np reaction at four selected energies. Dotted and solid curves are obtained 
in IA and the full model, respectively, with the MAID07 analysis and the amplitudes Ai. Shaded areas are our predictions with 
different parameterizations of the elementary 7 A — > irN amplitude as described in the text. Data are from Ref. [ll|]. Only 
statistical errors are shown. 



Results for the unpolarized differential yield p^|) for the 7c? — > ir°np reaction at four selected energies from 301 to 
455 MeV with a step of about 50 MeV are shown in Fig. [3] One can see a sizable effect from FSI, which is mainly due 
to np FSI. The yield is reproduced very well at the proton angles Q p > 40°. However, despite of the strong reduction 
of the unpolarized differential yield due to FSI, the full model clearly overestimates the data in the peak position. An 
analogous feature was found also in Ref. [16| where after integration over complete acceptance, i.e. with p™ m given 
by Eq. (fT2|) . we observed the predictions to lie well above data from Refs. 0,H3|- Similar overestimation, although in 
some variation of size, is present in all recent calculations (l7l l37l. l39j . 

It is difficult to find a reasonable explanation for this disagreement in the framework of the present model. First, 
reasons for a possible failure of the IA calculation are not apparent. The contribution of IA to the total cross section 
stems from the quasi- free ir°N process, which is evaluated very reliably. Second, as explained in Ref. [T^|, the NN FSI 
contribution comes mainly from the kinematic domains where nucleons N and N' in Fig. [TJb) are close to their mass 
shells so that the on-shell parameterization of the elementary operator T^n^ttN can be used in Eq. (f^Uj) . Because 
modern parameterizations of T^n^^n and T/vw are firmly established, we think that the effect of NN FSI is well 
under control. Therefore, the noticeable disagreement with the data cannot be explained by an underestimation of 
the JVJV FSI contribution. We also suppose that irN FSI is evaluated quite reasonably and it can not be responsible 
for the disagreement, especially if one takes into account that the wN FSI effect is very small. 

As discussed in Refs. [ill, H3, a possible way to resolve the problem is to take into account the absorption of the 
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produced pions by the pair of nucleons. Rough estimations of the absorption effect performed in Ref. [llj show that 
it can comprise more than 10% of the total cross section, thus leading to visible attenuation of the production rate. 

Model uncertainties of the predictions stemming from the choice of the elementary operator T-jv-nriV ar e shown in 
the shaded areas in Fig. [3] These areas have been generated, first, by results with the MAID07 and SAID (solution 
SP09K) multipole analyses and, second, with two sets of the invariant amplitudes, Ai and A' iy as described in the 
previous section. One can see that the model dependence is quite small, and even when taken into account, it does 
not resolve the problem with the description of the data at P < 40°. 
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FIG. 4: Unpolarized differential cross section for the 'yd — > n~pp reaction at four selected energies. Shaded areas are our 
predictions with different parameterizations of the elementary 77V" — > irN amplitude. Data are from Refs. [ll|] (•) and (°). 
Only statistical errors are shown. 



Predictions for the unpolarized differential cross section in the 7c? — > ir~pp reaction are compared in Fig. 0] to the 
data from Refs. [ill, El at the same energies as in Fig. [3] The agreement is quite reasonable, although the slight 
overestimation of the data is seen. The model dependence of the predictions is mainly due to the choice of the 
multipole analysis and is noticeable only at the lowest energy. It is compatible with the data errors of Ref. [4p| . 
Together with increasing E 1 this dependence diminishes and becomes quite small above 350 MeV. 

As is seen in Fig. [SJ the agreement of the predictions and the unpolarized ir + nn data is almost perfect. Again, the 
visible dependence of the results on the choice of the multipole analysis occurs only at 301 MeV. 

Taking into account the results for the unpolarized differential cross sections in the charged channels, one could 
expect the sum of the total cross sections for these two channels to be in reasonable agreement with the data but 
with a possible slight overestimation because of the above overestimation in the ir~pp reaction. However, as is seen 
in Fig. the deviation from the data of the GDH and A2 collaborations Q at E 1 > 300 MeV is more serious. At 
present, we are not aware of reasons for this disagreement. Figure [6] shows also that our results are in good agreement 
with these of AFS up to the peak position at about 320 MeV, but at higher energies our cross sections overestimate 
the AFS results. 
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FIG. 5: Same as Fig. 21 but for the yd — > ir + nn reaction. 
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FIG. 6: Unpolarized total cross section for the charged channels 7c! — > tt^NN. Shaded area is our predictions with different 
parameterizations of the elementary 7./V — > ttN amplitude. Solid curve is the AFS result 0. Data are from Ref. 0. 



B. Polarized yield and cross section 

Next we discuss the polarized yield and cross section. Figure [7] shows the helicity-dependent yield difference 

AY = dYp _ dYA (23) 
cj?^ - ^^ d/^l/jj ct^^i<p 

for the TT°np channel. Similar to the unpolarized yield, one can see a sizable reduction of AY/dfl p in IA due to FSI. 
However, the striking distinction from the unpolarized case is that now the model reasonably reproduces all the data 
not only those at O p > 40°. This might indicate that the above overestimation in the unpolarized yield is due to an 
overestimation in both dYp/dfl p and dYA/d£l p , but it disappears in the difference AY/dfl p . 
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FIG. 7: Helicity-dependent differential yield difference AY/dQ. p for the ■yd — > -K°np reaction. Dotted and solid curves are 
obtained in IA and the full model, respectively, with the amplitudes Ai and the MAID07 analysis. Shaded areas are our 
predictions with different parameterizations of the elementary 7iV — ► ttN amplitude. Data are from Ref. Only statistical 
errors are shown. 



It is seen in Fig. [7] that the dependence of AY/dVt p on the elementary operator is weak below 350 MeV, but it 
becomes quite visible at higher energies. 

Figures |8] and |9] show the helicity-dependent differential cross section difference 

Act _ dap _ da A ^ 

(l/^h'jY C^^~^yr- d/^h^ 

for the ir~pp and ir + nn channels, respectively. There is a good agreement with our results within the uncertainties of 
the theoretical predictions. Only a few ir~pp data points are underestimated at O w = 125°. The above uncertainties 
are quite significant only in the ir + nn channel and just in the angular region overlapped in the experiment 11]. 

An integration of Eq. ([2"4"]) over the solid angle W gives the difference Aa = ap — a a of the total cross sections. 
This difference for the semi-exclusive channels jd — > ^NN and jd — > tt°X has been measured in Refs. 0, [l(| at 
energies from 200 to 430 MeV. Because we do not have a model for the coherent reaction 7c? —> ir°d, we are able 
to make predictions for the second channel. Our results for the tt NN channel are shown in Fig. [TU] Within 
theoretical uncertainties, the model satisfactory describes the data of Ref. [9] below 300 MeV, but there exists 
noticeable overestimation at higher energies, i.e. one observes the situation analogous to that for the unpolarized 
cross section. Reasonable agreement with the data of Ref. [l(| is achieved except for two points in the peak region. 
Our results are in quite good agreement with the AFS predictions but overestimate these from Ref. [13] at E 1 > 280 
MeV. 

Having measured the unpolarized yield and cross section as well as AY/dil p and Act /dQ n , the authors of Ref. [ll[ 
have attempted to separate the parallel and antiparallel components of the yield and cross section. To do this, they 
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nN amplitude. Data are from Ref. [ill ]. Only statistical 



da + (25) 



used the follows relations: 

dn 3 v~ dn ' dn 

da _ V2fd(T\\ da Q \ 

where an and <jq are the components of the cross section corresponding to the deuteron spin states = ±1 and 
rrid = 0, respectively, and T20 is the tensor target asymmetry 0, Hl|. Note that the component an is the average of 
the parallel and antiparallel cross sections 

da\\ = l/da P da A \ 

dn 2 V dfl dSl J' [ ' 



It is follows from Eqs. and that 



da = V2 dgy 



Evaluations of T20 in Refs. [la . [l7| showed that in the angular regions accessed in conditions of the experimental 
setup of Ref. [ll|, i.e. at 0^ > 20°, the absolute value of T 2 o does not exceed 0.05. Therefore, one has the following 
approximate relations for 0^ > 20°: 
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FIG. 9: Same as Fig. [5] but for the 7c? —¥ -K + nn reaction. 
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FIG. 10: Difference Act for the charged channels 7c! — > tt^NN. Shaded area is our predictions with different parametrizations 
of the elementary 7 TV — > ttN amplitude. Solid and dashed curves are the AFS [3] and Schwamb [37| results, respectively. Data 
are from Refs. [9| (o) and [10|] (•). Only statistical errors are shown. 



Using Eqs. J24|), ([25]), (|27 |) . and (|29 |l one obtains 

dap 1 A _ da 1 Act 

dn ~ dn 2dn' K ' 

It should be noted that the semi-inclusive asymmetry T20 resulting from the integration of the corresponding right- 
hand sides over the kinematic regions appearing in Eqs. P^]) - ([To r ]) was not evaluated in Refs. 0,113 • However, we 
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have checked that such an integration gives IT20 1 < 0.08 at 25° < P < 65° and, therefore, the relations similar to 
Eqs. J2H]) and (|3T)|) are valid for the measured yields as well. 

Figures [TTHT31 demonstrate our predictions for the parallel and antiparallel yields and cross sections. A comparison 
with the experimental results from Ref. [TTJ] is also presented in these figures. One can see that the description of the 
data is quite reasonable, but not perfect. Such an agreement could be anticipated for the charged channels, but it is 
rather unexpected for the ir°np reaction. Currently, we are not aware of the reasons why the model that describes 
reasonably well the P and A components of the polarized yield, nevertheless, fails in reproducing the unpolarizcd 
yield. 
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FIG. 11: Helicity-dependent differential yields dYp/dQ p and dYA/dQ. p for the 7c! — > ir°np reaction. Solid and dashed curves 
are our predictions for dYp/dQ p and dYA/dQ, p , respectively, with the amplitudes At and the MAID07 analysis. Data are from 
Ref. pH for the parallel (•) and antiparallel (o) components. Only statistical errors are shown. 

We do not discuss here the shape of the angular behavior of the P and A components in the yield and differential 
cross section because the corresponding analysis has already been performed in Ref. [Tll |. As shown in that work, 
this behavior can be interpreted in terms of pion photoproduction on single nucleons, and it is a consequence of the 
dominance of quasi-free reaction mechanisms. 

C. The jd nNN contribution to the GDH integral for the deuteron 

Let us discuss now the contribution of the jd — > ttNN reaction to the deuteron GDH integral and compare our 
results with those obtained in Refs. [IH, [l5j|. Figure [T4l p resents the difference Act as a function of the photon energy 
Ej for three channels. As was shown in Refs. (l3l. Il5j. the behavior of the difference repeats mainly that for the 
elementary reaction 77V — > nN with some distortion due to nuclear effects. For instance, the negative value of Act in 
the charged channels at threshold energies is because of the dominance of the Eq + multipole producing the antiparallel 
transition. The M\ + multipole provides the strong positive value in the A region that Fig. 1141 clearly demonstrates. 
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FIG. 12: Helicity-dependent differential cross sections dap/dQ-n and daA/d^l^ for the 7<i — > n~pp reaction. Solid and dashed 
curves are our predictions for dap/dO.-^ and daA/dQ.^, respectively, with the amplitudes Ai and the MAID07 analysis. Data 
are from Ref. [Tj| for the parallel (•) and antiparallel (o) components. Only statistical errors are shown. 

As seen in Fig. [TJJ our results are quite close to those obtained in Ref. Visible difference may be seen only in 
the neutral channel in the first resonance region. 

The above disagreement also manifests itself in the value of the GDH integral, integrated up to 1.5 GeV. The 
upper limit of 1.5 GeV is chosen in accordance with an option from Ref. [l3[. Taking into account the theoretical 
uncertainties, we obtain the values of 241 ± 4, —20 ± 5, and 15 ± 20 (in units of /ib used also below) for the TT°np, 
ir~pp, and ir + nn reactions, respectively, to be compared to the AFS predictions, respectively, 222, 1 , —18.94, and 
2.51. One can see that our value for the w°np channel exceeds that of AFS by about 10%. Note here that the use 
of the SAID parametrization allows one to extend the upper limit in the GDH integral up to 2 GeV. With such an 
extension we have not found any visible changes in the result for 1.5 GeV. This means that the irNN term in the 
deuteron GDH integral essentially reaches convergence at these energies. 

The evaluation of the GDH integral, integrated up to 350 MeV, has been performed in Ref. [l5[ . Results of that 
work with the MAID03 model for the elementary jN —> irN operator are, respectively, 123, —29, and —9 for the 
7r°np, ir~pp, and ir + nn reactions, which are in complete agreement with our predictions of 123 ± 2, —34 ± 7, and 
— 7 ± 10, respectively. 

We have also studied the influence of the different parametrizations for the elementary nN amplitude on 

the value of the GDH integral, integrated to the maximum energy of 1.65 GeV. As is seen in Table HI the theoretical 
uncertainties in the ir°np reaction are quite small. They are more pronounced in the charged channels, especially in 
the Tr + nn one. This considerable dependence of the results is mainly due to the multipole analyses, but the effect of 



1 In fact, this value has not been given explicitly in Ref. |13| . We have extracted it from Fig. 5 of that work making use of a digitation 
procedure. 
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FIG. 13: Same as Fig. 1121 but for the yd — » ix + nn reaction. 
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FIG. 14: Difference Act for the separate channels of the yd — > nNN reaction. Shaded areas are our predictions with different 
parametrizations of the elementary yN — > irN amplitude. Curves are the AFS predictions 



a choice of the invariant amplitudes is also quite visible. Note that the noticeable dependence of the nucleon GDH 
integral evaluated with the SAID and MAID analyses was found in Refs. [l8l l42j. The sensitivity of the yd — > wNN 
contribution to the deuteron GDH integral, integrated up to 350 MeV, to a choice of the elementary operator have 
been studied in Ref. [l5j], and considerable dependence of the results obtained with the ELA and MAID03 operators 
have been found. 

Keeping in mind all the numbers given in Table HI we can quote our final value as 235 ± 25 fib for the contribution 
from the wNN channel to the GDH integral for the deuteron, integrated up to 1.65 GeV. As seen from Table HI the 
main part of the uncertainties stems from the ir + nn channel. It is worth mentioning that within the uncertainties of 
±25 /ib, the AFS result of 27.31 fjb [l3( for the total contribution to this integral for the deuteron is compatible with 



16 



TABLE I: Sensitivity of the GDH integral (in /tb), integrated up to 1.65 GeV in the 7c! — > ttNN reaction to the parametrization 
of the elementary 'yN —¥ nN amplitude. 



Parametrization 


Tv°np 


7T pp 


ir + nn 


tvNN 


MAID, Ai 


235 


-16 


6 


225 


MAID, 


238 


-24 


-4 


210 


SAID, Ai 


239 


-14 


35 


260 


SAID, A'i 


241 


-22 


25 


244 



the GDH sum rule prediction of 0.65 fjh. 

V. CONCLUSION 

The central aim of the present work is the study of the 7c? — > ttNN reaction with the polarized photon and deuteron 
and its contribution to the GDH sum rule for the deuteron. To evaluate the reaction amplitude, we have used the 
diagrammatic model elaborated previously in Ref. [161 ] . which takes into account diagrams corresponding to IA as well 
as NN and ttN interactions in the final state. For the elementary operator of pion photoproduction on the nucleon, 
its on-shell form generated by the recent multipole analyses, SAID and MAID, has been used. Particular emphasis 
has been placed on the discussion of possible uncertainties introduced into the model by this operator. 

Comparing our results with the recent experimental data obtained by the GDH and A2 collaborations in the energy 
domain from 300 to 450 MeV 0-[H|, we have rediscovered a known problem consisting in the sizable overestimation 
of the unpolarized data for the neutral channel. At the same time, the unpolarized data for the charged channels are 
reasonably reproduced. 

We have found that the helicity-dependent differential yield difference AY/ dfl p for the neutral channel and differ- 
ential cross section difference Aa/dfl^ for the charged channels are reproduced fairly well. Also, the model provides 
a quite satisfactory description of the separate parallel and antiparallel yields and cross sections in all the channels. 

We have also studied the uncertainties arising from the form of the elementary operator of pion production on the 
nucleon. Regarding the unpolarized observables, the sizable dependence has been found only for the charged channels 
at the lowest energy of about 300 MeV. The difference Ac/df^ manifests the weak sensitivity to the operator in 
the ir~pp channel in the full energy domain, and the difference AY/dtt p in the Tr°np channel does the same below 
350 MeV. However, a large sensitivity has been found in the ir + nn reaction from 300 to 450 MeV at 45° < 0^ < 125°. 

Within our model, we have evaluated the contribution of the irNN channel to the GDH integral for the deuteron, 
integrated up to 1.65 GeV. In addition, the influence of the elementary operator on its value has been estimated. Our 
final result is 235 ± 25 /ib. The uncertainty of ±25 /ib originates mainly from the ir + nn channel. Accepting these 
values as error bars for the value of 27.31 /ib obtained in the AFS framework 13] for the deuteron GDH integral, one 
can conclude that within these error bars, the AFS result is compatible with the GDH sum rule prediction of 0.65 /ib. 

Further improvements of the present model are needed to resolve, at last, the problem with the strong overestimation 
of the unpolarized cross section in the w°np channel. In particular, one should evaluate contributions of the two-loop 
diagrams which take into account effectively the absorption of the photoproduced n meson by the NN pair. On the 
experimental side, it would be very important to continue polarization measurements in the "fd — > TtNN channel. 
Specifically, data on the tensor asymmetry T20 together with these on the unpolarized cross section could provide the 
currently absent information on the rrid — component of the reaction amplitude. 
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